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Abstract. - We relate the Fermi-Dirac statistics of an ideal Fermi gas in a harmonic trap 
to partitions of given integers into distinct parts, studied in number theory. Using methods 
of quantum statistical physics we derive analytic expressions for cumulants of the probability 
distribution of the number of different partitions. 



Introduction. - Standard textbooks use grand canonical ensemble to describe ideal quan- 
tum gases [1]. It has been pointed out by Grossmann and Holthaus [2], and Wilkens et al. [3] 
that such an approach leads to pathological results for fluctuations of particle numbers in 
Bose- Einstein condensates (BEG). This has stimulated intensive studies of BEG fluctations 
using canonical and microcanical ensembles for ideal [4, 5] and weakly interacting [6] gases. 
Grossmann and Holthaus [2] related the problem for an ideal Bose gas in an harmonic trap 
to the number theoretical studies of partitions of an integer £ into M integers [7] , and to the 
famous Hardy- Ramanuj an formula [8]. In the series of beautiful papers, they have been able 
to apply methods of quantum statistics to derive non-trivial number theoretical results (see [9] 
and references therein). In this Letter we apply the approach of Ref. [9] to Fermi gases, and 
relate the problem of an ideal Fermi gas in an harmonic trap to studies of partitions of £ 
into M distinct integers. We show that the probability distribution of the number of different 
partitions is asymptotically Gaussian and we calculate analytically its cumulants. Our results 
are complementary to those of Tran [10], who has calculate particle number fluctuations in 
the ground state (Fermi sea) for the flxed energy and number of particles. 

Ideal Fermi gas. - We consider a microcanonical ensemble of N thermally isolated, 
spinless, non- interacting fermions, trapped in a harmonic potential with frequency loq. The 
discrete one-particle states have energies ty = hwf^iy + 1/2), = 0, 1, . . . . A microstate of the 
system is described by the sequence of the occupation numbers {n^}, where Un = 0,1 and 
riy = A. The total energy of a microstate {ny\ can be written as: 
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E{{n,}) = Ne^ + huj^£, 
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where integer: 

oo 

E^^n^v (2) 

determines the contribution to E from one-particle excited states. Note, that f is a sum 
of M distinct integers, where M = 7V— lifno = 1, orM = N otherwise. Hence, in 
order to construct the microcanonical partition function T{E,N), one has to calculate the 
number of partitions of the given integer £ into M distinct parts. This, in turn, is a standard 
problem in combinatorics [11]. If we call the number of distinct partitions <&(£,M), then 
T{E, N) = ^{£, M) and the physical problem of finding r{E, N) is mapped to a combinatorial 
one. 

Introduce the total number of distinct decompositions: 

A/™ ax 

£ HS,M), (3) 
and consider the probability (relative frequency): 

that exactly M distinct terms occur in the decomposition of £. This is the central object 
of our study. In the case of bosons Pmc{£, M) possesses clear physical interpretation, as it 
was explained in Ref. [9]. This is the probability of finding M excited particles, when the 
total excitation energy is HloqE. For fermions, however, the ground state of the whole system 
corresponds to the filled Fermi sea. Since the number of partitions M differs at most by 1 
from the number of particles TV, we may regard I0J as a probability density in the fictitious 
Maxwell demon ensemble [5] with vanishing chemical potential. The distribution Q describes 
at the same time an interesting mathematical problem, which, up to our knowledge, has never 
been treated using physical methods. _ 

Using standard results from combinatorics [11], the exact expression for <l'(£,Af) can be 
written down: 

$(£-,M)^p(^£- Q-^),Af), (5) 

in terms of a recursively defined function P: 

for fc = or k > n 
P{n,k) = { P{n~l,k-l)+P{n-k,k) . (6) 

1 for k ^ n 

The above recursion enables straightforward, although tedious, numerical studies of the dis- 
tribution ^ (see Fig. [TJ. 

Canonical distribution. - To proceed further with the analytical study oipmc{£TM) we 
introduce the corresponding canonical description. Instead of fixing the energy E of the each 
member of the ensemble, we consider an ensemble of systems in a thermal equilibrium with a 
heat bath of temperature T. The canonical partition function Z{(], N), (3 = 1/kT, generating 
the statistical description in this case, can be written using as: 

oo 

E £=0 
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P (e,M) 




Fig. 1 - The distribution p™c(f , M) for E = 1000, 2000, 3000, and 4500. 

where f = term vanishes. Introducing h — fkuoP we can define a combinatorial partition 
function: 



(8) 



£ = 1 



which can be viewed as the generating function of the microcanonical quantities $(£,M). 
The canonical probability, corresponding to is then given by: 



Pcn{b,M) 



Zcib,M) 



S=i' 



-b£ 



(9) 



A convenient and physically relevant way of an analytical description of the distribution © is 
provided by the cumulants. For a probability distribution {p{M), M > 0} of an integer- valued 
random variable, kth cumulant k''^'^ is defined as [12]: 



kC") = ( z— ] Inp(z) 



(10) 



where p{z) = J2m'^ p{M) is the (real) characteristic function of the distribution. The 
relation of the first few culumants to the central moments /x^'^' is the following: 

k(i)=^(i)^M k(2)=m(2) «:(3)^^(3) «:(4)=^(4)_3(Ai(2))2 , (11) 



where M = Sa/ ^^"^ mean value. Denote as 'E,c{h,z) = J2Ti=o ^c{b, M) = 
]^j^q(1 + ze^'"^) the grand partition function, corresponding to the canonical partition func- 
tion ijHl with z being an analog of fugacity. We are interested in excited part so we take not 
Sc(6, z) but rather Eex{b,z). 



Sea,(5,z) = n (1 + ze- 



(12) 
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We substitute ((HJ) in to the definition ifTUI) and obtain: 

(13) 



Now, we can apply the procedure of Ref. [9] in order to obtain a compact expression for 
the cumulants nin (b) as well as their asymptotic behavior for 6 <C 1. First, we note that: 

oo 

lnSe,(6,z) = ^ln(l + ze-^'^). (14) 
1^=1 

and we can always find such a neighborhood oi z — 1 that ze^^ < 1. Since, according to 
H13|l . we are only interested in the behavior of Sea: in the vicinity of z = 1, the logarithms in 
the second term can be Taylor expanded; using then the Mellin transformation of the Euler 
gamma function T{t) [13], we obtain the following integral representation of InSej, : 



1 



r+lOO 



In Se, (6, z) = — dt h-\{t)T{t)ft+,{z) , (15) 

ZTTl 



where C,(t) is the Riemann zeta function [12], and 

(-1 



/"W = >.^^^- (16) 



n=0 



is the fermionic function [1]. The parameter r > was chosen in such a way, that the 
integration line in (|15|l is to the right from the last pole of the integrand. 

Finally, we substitute Eq. (|15|l into (|13ll and with the help of the following properties of 
the fermionic function /„: z dfa{z)/dz ~ fa~i{z), /a(l) = (1 — 2"'^~")^(a), we find that: 

= ^ii,osonib) - 2'=4t!fco.oJ26) , (17) 

where 



'^S.o.o„(^) = ^ / dt b'^mcmt + 1 - fc) (18) 



r +100 



-loo 



are the canonical cumulants of the bosonic counterpart of the distribution @. 

Equation (|17|l is of central importance for this letter. It links the fermionic cumulants to 
that corresponding for bosons. The latter has been studied in detail in Ref. [9] . In particular [9] 
gives the asymptotic expressions of the first few bosonic cumulants for 6 <C 1. Using those 
together with H17I18|I allow us to find the asymptotic expressions for the cumulants of the 
fermionic canonical distribution ^ : K^^ib) = 7rVl2& + In 2/2 + 6/24 + ©(foi^-S), nil^ib) = 
ln2/&-l/4 + 6/48 + C'(&2-5), K'-^}{b) = l/2b- 1/8 + 0{b^^-^), nf^ib) = l/4fe- 6/96 + ©(fe^-'"^), 
i^mib) = —1/16 + 0{b^^-^), etc. The dependence on b in the fermionic case is of lower order 
than the correspponding one in the bosonic case. For example nin is constant up to the terms 
of the order b 



11.5 



Microcanonical distribution. - We come back to the study of the distribution Our 
aim now is to find the cumulants Kml{£) of this distribution. Applying the definition UlUf) to 
we find that: 



nl^li£)={zj-Y\nT{£,z) 



(19) 
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where: 



(20) 



M=0 



is another generating function of $(£, A/), complementary to Zc{b^ M). 
In order to calculate T(£, z), we first note, that: 

oc oo oo 

M=0 £=1 £=1 

Introducing a new variable x = and witting Sex(— Inx, z) = Sea;(a;, 0) we obtain: 



(21) 



(22) 



£=1 



that is, T(iS',z) are the coefficients in the power series expansion of ^6^(2;, z) w.r.t. x. We 
treat a; as a complex variable and use an integral identity: 



ZTTl / X'^ 



(23) 



where the integration contour is any closed loop surrounding the origin, oriented anti-clockwise. 
The integral in H2^^|l can be approximately evaluated using the saddle point method similar to 
that used in Refs. [1,9]. The implicit equation for the saddle point X{){z) is: 



d 

£ + 1 = -T^lnSea;(6o(2),z) , 
00 

with 60(2) = — In 2:0(2:). From Eq. (|24|) we obtain in the Gaussian approximation: 



(24) 



lnT(f,z) « ln.^e2.(6o(z),z) +£60(2) - iln27r- i In 



d 

db ' ^^-ex{ba{z),z) 



(25) 



The above solution for T(£, z) is of course similar to the one obtained in a bosonic case in 
Ref. [9], the difference being just in the form of the grand canonical partition function Se^. 
The microcanonical cumulants are then calculated from H25|l according to the equation 



lnT(£,z) 



(26) 



z=l 



where from Eq. dbo{z)/d2 



— — (5o(l))/9^Kin (60(1)), where in the last expres- 
sion we used the definitions of the canonical cumulants (|13|l . Generally, from (|13|l . (|25|l and 
H13|l it follows that the microcanonical cumulants Kml^S) can be expressed in terms of the 
canonical cumulants (6) and their derivatives w.r.t. the temperature parameter b. For 
example for the first cumulant Kmc one obtains: 



4'i(f)=4n^(&0(l)) 



152k(J,)(6o(1)) , dbK^^Jiboil)) , ia3K(°)(6o(l)) 



-I- 



1 + 



2 92^(0) (60(1)) d^K^^J (&o(l)) V 2 (feo(l)) 



(27) 
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For 6 <C 1 we can use calculate the first few microcanonical cumulants explicitly: 



36(ln2)= 



V3 



[ - 367r2 In 2 + tt* + 432(ln 2)^] VS 



|7r4- 547r2 + 864(ln2)3 



(28) 
(29) 

(30) 



3%/3 



— [An^ In 2 + 2160(ln 2 f - 21671^ (In 2)^ + 3tt*] VS 

t:^ + 31104(ln2)'* + 2771"* - 25927r2(ln 2)^ + 367r^ln2 + ©(f^*) .(31) 



TT 

1 r 1 



16 



Using (|28I31|) we obtain analytic expressions for the skewness and excess [12] of the distri- 
bution 7i(£) = KmliS) / (umliS))^ , l2i£) = Kml{£) / {umliS))"^ ] thesc parameters measure 
the deviation of the distribution Q from the Gaussian one. We obtain: 



71 (f) = 
72(f) = 



0.12894 



4Vf 

1.2001 



(32) 



and hence in the limit £ —f 00 both 71 and 72 vanish, implying that for large integers 
£ the distribution Q approaches Gaussian. This is a different situation from the bosonic 
case, where, as shown in Ref. [9], both skewness and excess are non-zero in the limit of large 
integers, so that bosonic microcanonical distribution does not approach Gaussian. In Fig|21we 
compare the numerically obtained dependence of the skewness and excess with the analytical 
predictions of Eq. (|32|l . 

Summarizing, we have studied relations between the Fcrmi-Dirac statistics for an ideal, 
harmonically trapped Fermi gas, and the theory of partitions of integers into distinct parts. 
Using methods of quantum statistical physics, we have described analytically properties of the 
probability distribution of the number of different partitions. We anticipate that the results 
concerning microcanonical ensemble for fermions can be used to characterize degeneracies of 
fermionic many body states, that play an essential role in the process of laser cooling of small 
samples of atoms in microtraps in the so called Lamb-Dicke limit [14]. 
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